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Ian M. Hoffman 

Wittenberg University, Springfield, OH 45501 
Accepted to ApJ 

ABSTRACT 

We present the first astronomical detection of the 14 NH3 (J, K) = (10, 6) line: nonthermal emission 
at several velocities in the Galactic star-forming region NGC 7538. Using the VLA we have imaged 
the (10,6) and (9,6) ammonia masers at several positions within NGC 7538 IRS 1. The individual 
sources have angular sizes < 0.1 arcseconds corresponding to brightness temperatures Tb > 10 6 K. We 
apply the pumping model of Brown & Cragg, confirming the conjecture that multiple ortho-ammonia 
masers can occur with the same value of K. The positions and velocities of the (10,6) and (9,6) masers 
are modeled as motion in a possible disk or torus and are discussed in the context of recent models of 
the region. 

Subject headings: HII regions — ISM: individual (NGC 7538) — ISM: molecules — Masers — Radio 
continuum: ISM — Radio lines: ISM 



1. INTRODUCTION 

In an ongoing series of papers (Hoffman & Kim 2011a, 
hereafter Paper I; Hoffman & Kim 2011b, hereafter 
Paper II) we have examined the ammonia masers in 
NGC 7538 with the hope of developing new constraints 
on the poorly understood phenomenon of nonthermal 
emission from nonmctastable (J > K) states. Although 
thermal emission from metastable (J = K) states is well 
developed into a method of thermometry of the emit- 
ting volumes {e.g., Maret et al. 2009), nonmetastable 
transitions have received comparatively little attention. 
Indeed, nonmetastable ammonia masers were discovered 
serendipitously in 1984 by Madden et al. (1986) after 
which only one emitting region (W51) was studied in 
depth (e.g., Wilson, Johnston, & Henkel 1990; Pratap 
et al. 1991) culminating in a pump model by Brown and 
Cragg (1991, hereafter BC91) that has never been tested. 
With the exception of the discovery of a new maser source 
by Walsh et al. (2007), there has been no study of non- 
metastable ammonia masers for 20 years. 

Nonmetastable ammonia masers are found in Galac- 
tic high-mass star-forming regions. Only five masering 
regions are known: NGC 7538, W51, DR21(OH), W49, 
and NGC 63341. Although a large-scale, dedicated sur- 
vey for nonmetastable ammonia masers has never been 
conducted, the HOPS survey by Walsh et al. (2011) in- 
cluded the frequencies of the (J,K) — (11,9) and (8,6) 
14 NH3 transitions. With a detection threshold of ss 2 Jy 
(suitable to detect only the strongest known ammonia 
masers), the HOPS survey resulted in the detection of 
only one new source of possible ammonia-maser emis- 
sion despite surveying more than 500 Galactic sources 
of water-maser emission. In contrast, the sensitive (~ 
100 mjy) survey of 17 star-forming regions conducted by 
Madden ct al. (1986) resulted in four detections. Thus, 
the prevalence and diagnostic utility of ammonia masers 
in the study of star formation remain open questions. 

Brown and Cragg (BC91) developed a pumping scheme 
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for ortho-ammonia (K = 3n;n — 0,1,2,...) in which 
there is a relationship among the states having the same 
azimuthal angular momentum K but different total an- 
gular momentum J. These "rungs" of value J within a 
"ladder" of value K are suggested to be excited simulta- 
neously, resulting in emission from multiple rungs for cer- 
tain physical conditions. Astronomically, only masering 
in isolated rungs - such as (6,3) without accompanying 
emission in (5,3), (7,3), or elsewhere in the K = 3 ladder 
- has ever been observed from sources. One exception 
may be NGC 63341 in which Walsh et al. (2007) report 
possible maser emission in (10,9) and (9,9) accompany- 
ing the known maser in (11,9). In the laboratory, Willey 
et al. (1995) have pumped an isolated (4,3) maser. The 
suggested interrelationship of the ladder rungs is impor- 
tant for two reasons: (1) many of the ammonia masers 
in existence may yet be undetected since the frequencies 
of most rungs' transitions lie outside the range of tradi- 
tional radio astronomy receivers, and (2) the conditions 
for simultaneous excitation of multiple transitions are far 
tighter than for one transition and, therefore, more valu- 
able as diagnostics of stellar environments. 

With the study described in this paper, we have set 
out to improve both the frequency coverage and angu- 
lar resolution of the study of nonmetastable ammonia 
masers. With the advent of the upgraded Karl G. Jan- 
sky Very Large Array (VLA) of the NRACQ, receivers 
are now available at the frequencies of many transitions 
of interest. Furthermore, improving upon the angular 
resolution of single dishes with interferometric imaging 
provides three significant benefits to understanding the 
excitation mechanism: (1) the maser positions may be 
determined relative to the continuum emission in the re- 
gion in order to constrain the pumping environs, (2) the 
angular sizes of the masering regions may be constrained, 
providing for tighter limits on the brightness temperature 
of the emission and on the amplification gains, and (3) 

The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 1.— VLA spectra for (top) 14 NH 3 (10,6) and (bottom) (9,6) emission from NGC 7538 IRS 1. Both data sets have had the 
contribution from continuum emission subtracted in the u, v plane. For display in this figure, each of these data sets has been convolved 
with the beam from Paper I: 4"6 X 1'.'9 at a position angle of 34°. The gray vertical arrows in the (10,6) spectrum denote the velocity 
features whose sky positions are plotted in Figure [2] The gray curves in the (9,6) spectrum represent the fitted Gaussian functions that 
are summarized in Table [T] and whose sky positions are plotted in Figure \3\ 



sufficiently precise angular positions of different transi- 
tions and velocity components will aid in discriminating 
isolated from coincident pumping. In addition, the gen- 
eral benefit of any maser study is that the kinematics of 
the gas in the region can be determined at the scale of 
~ 100 AU. With a deeper understanding of these basic 
parameters, nonmetastable ammonia masers can be built 
into a robust tool for understanding star formation. 

2. OBSERVATIONS AND RESULTS 
2.1. Ku-band VLA 'C Configuration 

On 2012 January 31 we observed the 14 NH 3 (10,6) am- 
monia transition at i/ rcst = 16.31932 GHz, employing the 
15 VLA antennas retrofitted with new fTu-band receivers 
at the time. The relative separations between the us- 
able antennas ranged from 0.1 to 2.6 km, resulting in a 
circular synthesized beam of FWHM diameter 1.8 arc- 
seconds. The flux, phase, and bandpass calibrators used 
were 3C48, 2322+509, and 3C147 with approximately 17 
minutes of integration on the IRS 1 target. We recorded a 
single circular polarization over a 1-MHz bandwidth cen- 
tered in frequency on i>lsr — — 53.8 kms™ 1 divided into 
256 spectral channels resulting in a velocity resolution of 
0.1 kms -1 and a total velocity coverage of approximately 
18 kms" 1 . The data were reduced using AIPSQThc rms 
background noise in a channel image is 17mJy beam -1 , 

2 The Astronomical Image Processing System is documented at 
http : / /www . aips . nrao . edu/ 



consistent with instrumental expectations. 

In the top panel of Figure [T] is displayed the ammonia 
spectrum. Due to spatial and spectral blending, Gaus- 
sian components could not be fit reliably to the velocity 
profile. Nevertheless, obvious velocity features (denoted 
by vertical arrows in Figure [T]) were chosen subjectively 
in order to select the channel images from which the po- 
sitions in Figure [2] were determined. None of the velocity 
components are coincident on the sky, indicating that the 
emission arises in several distinct kinematic volumes. 

Most of the masers are not angularly resolved in 
the VLA observations and have an upper limit on de- 
convolved size of approximately 1 arcsecond. These 
sizes correspond to a lower limit on brightness tem- 
perature of Tb > 2000 K. The maser component at 
w lsr = ~ 59.9kms~ 1 has a resolved size of 40 milliarc- 
seconds corresponding to a brightness temperature of 
Tb = 1.1 x 10 6 K. Assuming that this maser amplifies the 
16.3-GHz background continuum emission from IRS 1 
following Tb — Tb g e~ T , we find the gain to be r = —4.3. 
The measured parameters for all of the (10,6) emission 
is summarized in Table [T] 

2.2. K-band VLA l B' Configuration 

On 2012 June 15 we observed the 14 NH3 (9,6) ammo- 
nia transition at v lcst — 18.49939 GHz, employing all 27 
antennas of the VLA. The relative separations between 
the antennas ranged from 0.3 to 11.1 km, resulting in a 
synthesized beam of / .'33 x 0'.'26 at a position angle of 
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Fig. 2.— Positions of the 14 NH3 (10,6) masers in NGC 7538 
IRS 1 from the VLA 'C configuration observations of 2012 Jan- 
uary 31. The synthesized beam is circular with a FWHM diameter 
of 1.8 arcseconds. The diagonal labels are the velocities of the 
features denoted in Figure [T] and summarized in Table [T] For the 
location of the maser in each imaged spectral channel, the error 
bars are the position uncertainty of a fitted two-dimensional Gaus- 
sian model. The filled hexagon is the location of the peak of the 
16.3-GHz continuum emission from which the po sitio n offsets in 
this data set are measured. As discussed in Section l2.3l these error 
bars are appropriate for determining the relative separation of the 
features but do not represent the precision of the absolute position 
registration. 

— 15°; the array was not sensitive to image features larger 
than approximately 10 arcseconds. The flux, phase, and 
bandpass calibrators used were 3C48, 2322+509, and 
3C147 with approximately 45 minutes of integration on 
the IRS 1 target. We recorded a single circular polar- 
ization and a 2-MHz bandwidth centered in frequency 
on i»lsr = — 53.8kms _1 divided into 256 spectral chan- 
nels resulting in a velocity resolution of approximately 
0.2 kms -1 and a total velocity coverage of approximately 
32kms~ 1 . The data were reduced using AIPS. The rms 
background noise in a channel image is 14mJybeam _1 , 
consistent with instrumental expectations. 

Using the VLA, we detect the same velocity compo- 
nents as described in Paper II for the observations us- 
ing the Robert C. Byrd Green Bank Telescope (GBT), 
with the exception of Ulsr — ~~ 46.7kms _1 which is 
below the detection threshold of the current observa- 
tions. In the bottom panel of Figure [T] is displayed 
the VLA ammonia spectrum and fitted velocity compo- 
nents, which are summarized in Table [TJ Due to spa- 
tial and spectral blending, the fits are not reliable to 
the widths and amplitudes of the components having 
— 57kms _1 > wlsr > — SSkms" 1 . Nevertheless, the 
center velocities of the features are apparent and were 
used to select the channel images from which the posi- 
tions in Figure [3] were determined. 

The current results are also consistent with the VLA 
observations presented in Paper I, with the maser emis- 
sion near ulsr = — 60.2kms _1 lying southeast of the 
peak of continuum emission and of the masers near 
u lsr = — 57kms . Furthermore, all of the velocity 
components have positions at the location of IRS 1 in 
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Fig. 3.— Positions of the 14 NH 3 (9,6) masers in NGC 7538 
IRS 1 from the VLA 'B' configuration observations of 2012 June 
15. The images have a synthesized beam of 0"33 X 0''26 at a 
position angle of —15°. The diagonal labels are the velocities of 
the features summarized in Table [l] and plotted in Figure [l] For 
the location of the maser in each imaged spectral channel, the 
error bars are the position uncertainty of a fitted two-dimensional 
Gaussian model. The filled diamond is the location of the peak of 
the 18.5-GHz continuum emission from which the posit ion offsets 
in this data set are measured. As discussed in Section 12.31 these 
error bars are appropriate for determining the relative separation 
of the features but do not represent the precision of the absolute 
position registration. The gray dashed line de notes the direction 
of the Keplerian gradient discussed in Section 13.21 and displayed 
in Figure [4] The star symbol is the fitted location of the central 
object of the Keplerian model, with error bars denoting formal 
fitting errors. 



NGC 7538. These results are consistent with the GBT 
observations presented in Paper II and allow us to rule 
out the association of these ammonia masers with any 
other sources in the region (such as IRS 9 or MM2). 

The velocity feature near i>lsr = — 57kms~ 1 is far 
less pronounced in the current spectrum than in previous 
observations, indicating some variability among the sev- 
eral blended components between — 57kms~ 1 > ulsr > 
— 53kms . In addition, the entire amplitude scale is 
w 2 times larger in the current paper than in Papers I 
and II. The flux scale of the current observations was 
confirmed independently against both 3C48 and 3C147 
with consistent results, leaving the discrepancy an open 
question; we hesitate to conclude that all of the maser 
components varied in unison by a factor of two since mul- 
tiplicative scaling is a common instrumental uncertainty. 

As summarized in Table [U many of the masers are an- 
gularly resolved, having sizes on the order of 100 milliarc- 
seconds (« 300 AU at a distance of 2.65 kpc, Moscadelli 
et al. 2009). These brightness temperatures of Tg ~ 
10 6 K correspond to unsaturated amplification gains of 
t w —5 of the background continuum emission. 

As reported in Papers I and II, the velocity components 
near fLSR = — 60kms _1 have been relatively unchanged 
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TABLE 1 

Fitted Properties of Ammonia Emission 



"LSR 


/ 


AtlpwHM 


^mascr 


T B 




T 


(kms" 1 ) 


(Jybm- 1 ) 


(kms -1 ) 


(mas) 


(10 6 K) 






14 NH 3 (10,6) 


-61.4 a 


b 


b 


<1000 


>0.001 


< 


— 1.1 


-59.95(1) 


0.37(1) 


0.55(3) 


40 


1.1 




— 4.3 


-56. 7 a 


b 


b 


<800 


>0.002 


< 


-1.7 


-55. 8 a 


b 


b 


<900 


>0.003 


< 


-2.0 


-54.9 a 


b 


b 


<600 


>0.005 


< 


-2.8 


-53. 9 a 


b 


b 


<700 


>0.002 


< 


-1.9 


i4 NH 3 (9,6) 


-60.20(4) 


1.78(8) 


0.63(6) 


70 


1.3 




-5.0 


-59.77(3) 


1.3(3) 


0.45(3) 


<90 


>0.6 


< 


-4.3 


-56.7(5) 


b 


b 


50 


2.0 b 




-5.5 b 


-56.1(2) 


b 


b 


60 


0.7 b 




-4.4 b 


-55.4(8) 


b 


b 


60 


0.6 b 




-4.3 b 


-54.8(7) 


b 


b 


<90 


>0.07 b 


< 


-2.2 b 


-54.3(9) 


b 


b 


100 


0.13 b 




-2.9 b 


-50.07(3) 


0.22(2) 


0.54(6) 


<150 


>0.04 


< 


-1.9 


-48.64(5) 


0.11(2) 


0.5(1) 


<200 


>0.01 


< 


-0.7 



Note. — The number in parentheses is the uncertainty in the 
final digit. The fitted angular sizes of the masers have a precision 
of one significant digit, from which the brightness temperatures and 
gains are calculated. 
a Approximate velocity of feature. 
b Fit uncertain due to blending of features. 

since their discovery in 1984 (there are unpublished re- 
ports of variability by a factor of 2.5 in the year following 
discovery, see Schilke et al. 1991). In contrast to these 
older masers, all of the redward emission that has arisen 
subsequently is found to lie along the same northeast- 
southwest line (see Fig. This velocity distribution is 
discussed in Section [3T2l 

2.3. Comparison Between Observations 

The absolute position registration between the 16.3- 
GHz and 18.5-GHz observations is limited by interferom- 
eter phase calibrations to a precision of approximately 
200 milliarcseconds. Both observations contain an im- 
age of the continuum radiation from IRS 1: the fitted 
peak of the image of the 16.3-GHz continuum emission 
lies northwest (position angle —61°) of the fitted peak 
of the image of the 18.5-GHz continuum emission by 50 
milliarcseconds. We assume these peaks to be coinci- 
dent in position and use them to register the two ob- 
servations to each other: the maser positions plotted in 
Figures [5] and [3] are relative to the peak of the contin- 
uum emission at that frequency. The relative location 
of the masers to each other in a single observation is far 
more precise (« 5 mas) than an unregistered compari- 
son between observations. The location of the peak of 
the continuum emission in the 18.5-GHz observations is 
aj2ooo = 23 h 13 m 45?364, 6 ]20ao = 61° 28' 10734. 

Sandcll ct al. (2009) discuss extensively the morphol- 
ogy of the continuum emission from IRS 1. Based on 
observations with very high angular resolution, they de- 
scribe differences in morphology at the 150-mas scale be- 
tween 14.9 and 22.4 GHz. Given the poorer angular reso- 
lution and smaller difference in frequencies of our current 
observations, we consider our assumed 50-mas registra- 
tion to be sufficiently precise for our discussion. 

Alternatively, we may assume coincident positions for 



(10,6) and (9,6) ammonia masers with comparable veloc- 
ities (similar to the registration method used by Surcis et 
al. [2011] for water and methanol masers). Most notably, 
the brightest spectral features at Ulsr ~ —59.9 kms -1 
and i>lsr ~ — 56.7kms _1 have comparable relative sep- 
arations at each frequency. The possibility of spatially 
coincident maser emission from both (10,6) and (9,6) is 
addressed in the next section. 

3. DISCUSSION 
3.1. Pump Model 

Brown and Cragg (BC91) suggest a pumping scheme 
for the (6,3) maser in W51 that accounts for the lack of 
maser emission in the (5,3) and (4,3) rungs of the K = 3 
rotational ladder. The model requires radiative transi- 
tions from the vibrational ground state that are driven 
by an adjacent infrared source; vibrationally excited am- 
monia has been detected in NGC 7538 (Schilke et al. 
1990). The subsequent rapid radiative decays back into 
the vibrational ground state would invert all rungs of 
the ladder were the populations not mitigated by com- 
peting collisions which thermalize some of the transi- 
tions. For certain ranges of infrared excitation, ammo- 
nia column, collision frequency, and collision selectivity, 
an ortho-ammonia maser can be produced in a specific 
rung. Although the pumping scheme is heavily depen- 
dent on the cross-sections for selective collisions between 
ammonia and molecular hydrogen, the BC91 model is 
found to be insensitive to the input of either the rates 
extrapolated from helium calculations or the rates calcu- 
lated specifically for molecular hydrogen. Furthermore, 
BC91 note that the conditions for masers with K = 6 are 
much different than for K = 3, and so these kinematic 
volumes in IRS 1 should be devoid of K = 3 emission. 
Indeed, the (6,3) and (4,3) transitions have been studied 
in NGC 7538 and neither show maser emission (Madden 
et al. 1986; Schilke et al. 1991). 

In the laboratory, Willey et al. (1995) successfully 
generated a nonmetastable ortho-ammonia maser in the 
(4,3) transition (in isolation from [5,3] or [3,3] emission) 
using exclusively collisional pumping (see also, Schilke 
et al. 1991). In contrast to the analysis by Brown and 
Cragg, the laboratory maser was generated only in colli- 
sions with helium and could not be produced under any 
conditions with molecular hydrogen. Flower, Offer, and 
Schilke (1990) note that population inversion is expected 
to depend on the ratio of ortho-H2 to para-H2 partic- 
ipating in the collisions. Willey et al. (2000) went on 
to quantify significant differences between NH3-He col- 
lisions and NH3-H2 collisions. Furthermore, the labora- 
tory maser is explained using transitions from the K = 
ladder even though that process is ruled out of the BC91 
analysis. Nevertheless, all groups of authors recognize 
that the role of collisions in pumping is predicated upon 
having level- and parity-dependent selectivity. 

As noted by BC91, all theoretical work on ammonia 
masers has been limited by the lack of astronomical stud- 
ies of the rungs higher in energy than the masering rung. 
In the current paper, we present the first detection of 
masers in adjacent rungs of an ortho-ammonia ladder 
and discuss here briefly the implications of simultaneous 
maser emission from (10,6) and (9,6). We are currently 
undertaking observations toward a complete survey of 
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the ortho-ammonia ladders in NGC 7538 and withhold 
a larger discussion of level populations for a subsequent 
paper. 

Regardless of how the 16.3- and 18. 5- GHz observa- 
tions are registered (see §2.3j) . there is a position co- 
incidence of (10,6) and (9,6) maser emission at either 
u lsr ~ — 54kms~ 1 or at — 60kms _1 and — 57kms" 1 . 
Brown and Cragg explicitly address the possibility of the 
same volume of gas supporting both population inver- 
sions simultaneously. In the same way that they predict 
a (6,3) maser to be accompanied by (7,3) - and possibly 
also higher values of J - we have found the (9,6) maser 
to be accompanied by (10,6). Qualitatively, we confirm 
their conjecture that masers in a given ladder may not 
appear in isolation but rather in all rungs above a critical 
value of J = 9. Quantitatively, little can be concluded 
without the detailed cross-section information, especially 
the mitigation of the radiative decays that is necessary 
to produce these (10,6) and (9,6) results. With future 
work, a sufficiently complete census of the rungs in the 
K = Q ladder will allow the most selective level- and 
parity-dependent NH 3 -H 2 cross-sections to be deduced, 
even in the absence of laboratory results. 

3.2. Velocity Gradient 

In terms of both ages and positions, the (9,6) masers 
with wlsr ~ —60 kms" 1 are qualitatively different than 
the younger, redward masers having vlsr > —58 kms" 1 . 
Indeed, the straight-line position arrangement of the 
redward masers is accompanied by a monotonic veloc- 
ity gradient that is not shared by the older masers. 
The northeast-southwest orientation of the gradient (see 
Fig. \S§ is consistent with the orientation of many other 
proposed disks in the region (see Beuther et al. [2012] for 
a recent review). 

In our case, the gradient is well fit as Keplerian motion 
about an unseen source (the star symbol in Fig. [3]) that is 
significantly offset from the position of peak continuum 
emission. Figure 2] shows the fit to the gradient that re- 
quires a massive central object (180 ± 40M©) with a sys- 
temic velocity of vlsr ~ —70 kms" 1 . The masers occur 
only in the southwest half of the disk where the rotat- 
ing motion is away from the observer (counterclockwise 
as viewed from the northwest). This sense of rotation 
is opposite the other proposed disks with a northeast- 
southwest plane. The velocity of the fitted central ob- 
ject does not agree with the consensus velocity for IRS 1 
{vlsr ~ — 57 kms" 1 ). However, the fitted systemic ve- 
locity in our Keplerian model is comparable to that of 
the southern group of water masers described by Surcis 
et al. (2011) that lie approximately one arcsecond away 
from the fitted location shown in Figure [3J 

Although the fitted mass of the central object is unrea- 
sonably large in the Keplerian model, we hesitate to rule 
out this interpretation without understanding fully the 
contribution of the most redshifted masers to the model. 
Indeed, the most redshifted maser described in Paper II 
is not detected in the current observations. Since the 
highest velocities have the most leverage on the fit, we 
cannot conclude that a Keplerian model is unphysical 
without improved positions for all of the masers having 
u lsr > —51 kms" 1 . Nevertheless, this Keplerian model 
is only appropriate for the (9,6) masers, not the (10,6) 
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Fig. 4. — One of the two interpretations of th e velocity gradi- 
ent: the Keplerian model discussed in Section [3. 2 1 The data points 
are the masers whose positions lie along the gray line in Figure [3] 
with error bars denoting uncertainty in position along that line. 
The separations are measured from the fitted location of the cen- 
tral gravitating object, denoted with the star symbol in Figure [3] 
The conversion of angular position to the plotted linear distance 
assumes that the masers all lie with the central object in the plane 
of the sky at a distance of 2.65 kpc. The dashed line is the fitted 
gradient model. 

masers. 

As an alternative to a Keplerian interpretation, we 
consider that the (9,6) and (10,6) masers may mark the 
boundary of the torus/outflow region proposed by Surcis 
et al. (2011). Those authors model the velocity gradi- 
ents of methanol masers as arising in a torus with a flat 
rotation curve. The positions and velocities of the red- 
ward (9,6) ammonia masers agree well with the methanol 
gradients (their Figures 4 and 9). Furthermore, the dis- 
tribution of (10,6) masers (Fig. [2]) can also be understood 
in this manner. Using the parameters from Surcis et al. 
(see references therein) for the orientation of the torus 
(an outflow axis directed 40° west of north and an in- 
clination of 32°), the azimuthal position angles in the 
torus can be computed for the ammonia masers. We 
determine the position angles with the same sense and 
fiducial as presented in their Figure 9, and produce our 
analogous Figure [5] for ammonia. In fitting our data (the 
solid and dashed curves in Fig. [3]), we find a systemic 
velocity of wlsr = —55.4 ± 0.6 kms -1 and a rotation 
speed of 5.7 ± 1.1 kms -1 , in good agreement with the 
quoted methanol parameters of —57 kms" 1 and 7 kms" 1 , 
respectively. As discussed in Paper II, ammonia exhibits 
a higher range of velocities in IRS 1 than any other maser 
species besides water, and the ammonia masers at the 
ends of the velocity range are critical to our fit to the 
torus model. Although extraordinary for masers in the 
region, these ammonia velocities are nevertheless consis- 
tent with the range —65 kms" 1 < wlsr < —45 kms -1 
quoted by Klaassen et al. (2009) for the bulk rotation of 
SO2 gas in the torus. 

It is important to note that, in the torus model, all of 
the ammonia masers and some (group A) of the methanol 
masers lie in the part of the torus that is on the far side 
of the central object. This scenario requires that the 
ionized outflow in the foreground at those locations be 
optically thin at frequencies as low as 6 GHz (which is not 
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Fig. 5. — One of the two interpretatio ns o f the velocity gradi- 
ent: the torus model discussed in Section I3.2I The position angle 
is the azimuthal location of the masers in a suggested torus sur- 
rounding the central peak of continuum emission. The circles rep- 
resent (10,6) masers and the triangles represent (9,6) masers. The 
two gray triangle markers with PAfa —160° are the "older" (9,6) 
masers that have been excluded from the model. The horizontal 
error bars denote the uncertainty in the computed position angle 
due to the uncertainty in the sky positions of the masers. The 
solid curve is the best-fit torus model, with a systemic velocity of 
^LSR = — 55.4kms — 1 and a rotation speed of 5.7±1.1 kms" 1 , with 
the dashed curves representing the upper and l ower bounds of the 
fitted rotation speed. As discussed in Section 13.21 this modeling 
follows the analysis of Surcis et al. (2011) in the development of 
their Figure 9. 

inconsistent with observations, e.g. Sandell et al. 2009). 
Furthermore, the assumptions in Sections 12.11 and 12.21 
that the masers amplify the free-free continuum emission 
from IRS 1 are not valid and the maser gains are much 
stronger (more negative) than quoted in Table [TJ Also, 
note that the (9,6) masers near — 60 kms -1 have been 
excluded from the fit (the gray markers in Fig. [5J. 

One possible experiment to discriminate between a Ke- 
plerian interpretation and a torus interpretation of the 
ammonia gradient is to monitor the proper motion of the 
masers. In the case of Keplerian motion, the innermost 
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masers (the most redshifted) will complete an entire orbit 
in ps 50 years, and therefore vary on few-year timescales 
by either moving across the sky or vanishing as the beam- 
ing and velocity coherence become unfavorable. In the 
case of slow rotation with the torus, the masers would 
take approximately 200 years to move 100 milliarcsec- 
onds. 

4. CONCLUSION 

We present the first astronomical detection of the 
14 NH3 {J,K) = (10,6) line. Using interferometric im- 
ages of both the (10,6) and (9,6) masers in NGC 7538 
IRS 1, we find maser amplification gains of r w —5 and 
brightness temperatures of Tb ~ 10 6 K. Although many 
of the masers are isolated in both position and velocity, 
we note some possible coincidences between (10,6) and 
(9,6) emission which are discussed in the context of the 
pumping model of Brown and Cragg (1991). All of the 
new (9,6) masers reported in Papers I and II are found 
to have a velocity gradient. This gradient is consistent 
with both (1) Keplerian orbits about a very massive, un- 
seen object and (2) rotation in a torus around the known 
peak of continuum emission. The torus model can also 
be used to explain the positions of the (10,6) masers. In 
either case, the ammonia masers are apparently valuable 
markers of the kinematics of the molecular environment. 
Furthermore, the long-lived (9,6) masers originally dis- 
covered by Madden et al. (1986) do not fit well into either 
kinematic model and we leave open the possibility that 
they represent a different class than the other masers in 
this study. In an ongoing effort to quantify the value 
of nonmetastable ammonia masers in the study of star 
formation, we are undertaking additional observations of 
other rungs in the rotational ladders with both single- 
dish and interferometric observatories. 
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